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Abstract: In this paper, we analytically investigate multi-component Cahn-Hilliard and Allen- 
Cahn systems which are coupled with elasticity and uni-directional damage processes. The free 
energy of the system is of the form ^TVc : Vc -I- -I- c, z) dx with a 

polynomial or logarithmic chemical energy density an inhomogeneous elastic energy density 
IT®' and a quadratic structure of the gradient of the damage variable z. For the corresponding 
elastic Cahn-Hilliard and Allen-Cahn systems coupled with uni-directional damage processes, 
we present an appropriate notion of weak solutions and prove existence results based on certain 
regularization methods and a higher integrability result for the strain e. 
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1 Introduction 

Phase separation and damage are common phenomena in many fields, including material sciences, 
biology and chemical reactions. Such microstructural processes take place to reduce the total 
free energy, which may include bulk chemical energy, interfacial energy and elastic strain energy. 

The knowledge of the mechanisms inducing phase separation and damage processes is very im¬ 
portant for technological applications as for instance in the area of micro-electronics due to 
the ongoing miniaturization. The materials used in this area are typically alloys consisting of 
mixtures of several components (cf. |HCW9l] b 

Phase separation and damage processes are usually described by two separate models in the 
mathematical literature. To describe phase separation processes for alloys, phase-field models 
of Cahn-Hilliard and Allen-Cahn type coupled with elasticity are well adapted. On the other 
hand, damage processes for standard materials are often modelled as unilateral processes within 
a gradient-theory [FN96] . A phase-field approach which describes both phase separation and 
damage processes in a unifying model has been recently introduced in |HK11) . 

The main objective of this work is to prove under general assumptions existence results for 
multi-component systems where Cahn-Hilliard as well as Allen-Cahn equations are coupled with 
rate-dependent damage differential inclusions for elastic materials. We are interested in free en¬ 
ergies of the system which may contain a chemical energy of logarithmic or polynomial type, an 
inhomogeneous elastic energy and a quadratic term of the gradient of the damage variable. To 
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this end, we establish some regularization methods which enable us to show existence results for 
gradient terms jVzl^ of the damage variable 2; in the free energy even if the assumption p > n 
(n space dimension) is dropped. In contrast to |MR06[ IHKll) now the physical meaningful term 
|Vzp can be treated, cf. |FN961 lFre02) . In addition, we also provide a higher integrability result 
for the strain tensor. As a consequence, the chemical free energy may also have a logarithmic 
structure such that we are not restricted to polynomial growth as in [HKllj . We focus on the 
modelling of rate-dependent damage processes but we would like to mention that the results of 
this work can be extended to rate-independent systems (i.e. the dissipation potential is homo¬ 
geneous of degree one) by some modifications. In the following, we will introduce the model 
formally. 

The elastic material we want to consider in this work is an N component alloy occupying a 
bounded Lipschitz domain fl C M". To account for phase separation, deformation and damage 
processes in one model, a state at a fixed time point is described by a triple (u, c, z) where 
u : ft ^ K" denotes the deformation, c : —>■ the vector of the chemical concentrations and 

z : n —> K the damage variable. The meaning of the variables and its governing evolutionary 
process is explained more explicitly below. 

The mixture of the alloy is described by a phase field vector c = (ci,..., c^v), where the element 
Cfc for k = 1,... ,iV, denotes the concentration of the component k. Therefore, we will restrict 
the state space for c to the physically meaningful condition = I in 17. The constraint 

Cfc >0, k = 1,..., N, in 17 is also used for logarithmic chemical potentials (see below). 

If an alloy is cooled down below a critical temperature then usually spinodal decomposition and 
coarsening phenomena occur. Well established models for describing such effects are the Cahn- 
Hilliard and Allen-Cahn equations, which describe mass preserving and mass non-preserving 
phase separation in solids, cf. |Cah611 ILC821 IGur961 IAC791ICNC94) for model ling aspects. An¬ 
alytical investigations of Cahn-Hilliard equations can be found in [BCD^O^ IGarOOl ICMPOOl 
IGarOSal iGarOSbl IBP05) and for Allen-Cahn equations in |BdS04[ [CP081 ICGPGSlOj . respective¬ 
ly. The essential difference between these two equations is that the Cahn-Hilliard equation is 
a fourth order parabolic evolutionary equation expressible as a H~^ gradient flow of the free 
energy with respect to c whereas the Allen-Cahn equation is a second order parabolic equation 
arising from an gradient flow. More precisely, 

Allen-Cahn: dtc = -Mf - div(rVc) -b Wf{c) + IHe (e(w), c, z )), 

^ ' ' ' ( 1 ) 

Cahn-Hilliard: dtC = div(Mv( - div(rVc) -b Wf{c) + c, z))) . 

Here, denotes the chemical energy density, VP®* the elastic energy density, M the mobility 
matrix satisfying ^^^=1 = 0 for all fc = I,..., A^ and P the gradient energy tensor which is a 

fourth order symmetric and positive definite tensor, mapping matrices from into itself. 

In this work, VP®*' may be a chemical energy density of polynomial type, i.e. VP®*'(c) = VP®*'’P°*(c), 
or of logarithmic type, i.e. IP®*'(c) = VP®*'’*°s(c) (see (IA8I) 1. Note that phase separation only 
arises if the matrix A in (IA8I) is non-positive definite since the first term in (IA8II is convex. 

Elastic behaviour is modelled by a deformation variable u so that each material point x G ft 
from the reference configuration is located at x -b u{x). We use the assumption that the strain 
e is sufficiently small so that we can work with the linearized strain tensor given by e(u) = 
^(Vm -b (Vtt)‘). In this work, we will neglect inertia effects pil and volume forces 1. Therefore, 
the momentum balance equation div(cr) + I = pil from the continuum mechanics becomes a 
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quasi-static force equation, i.e. 


div((T) = 0. (2) 

The stress tensor a is defined by i.e. as the derivative of the elastic energy with respect to 
the strain. 

Analytical results for multi-component Cahn-Hilliard equations coupled with elastic deformations 
can be found in |Gar00] while Allen-Calm systems with elasticity are, for instance, studied in 
[BW05) . Finite element error estimates of Cahn-Hilliard equations with logarithmic free energies 
and concentration dependent mobilities are derived in |BB99) . Recent numerical results for 
Cahn-Hilliard and Allen-Cahn equations can be found in |BM10) . It turns out that different 
elastic moduli of the phases in the mixture influence the rate of coarsening and the morphology 
of the phases decisively [DM00] . Numerical investigations of elastic Cahn-Hilliard systems are 
conducted in [CRWOlj . 

The damage process we want to consider in this paper is uni-directional, i.e. it can only increase 
in time and the material is not able to heal itself. The phase field variable z satisfying 0 < 2 < 1 
is interpreted as damage in a way that z(x) = 1 stands for a non-damaged and z(x') = 0 for 
a maximally damaged material point a: G H. We assume that the damage in our model is not 
complete which means that a maximal damaged part has still elastic properties. These constraints 
lead to a differential inclusion formulation for the evolution of z which relates the derivative of 
the energy dissipation of the system depending on the rate of damage with the derivative of 
the free energy with respect to z. More precisely, we consider the doubly nonlinear differential 
inclusion (cf. |FN96j 'l 

0 e dp{dtz) - Az + W'^l{e(u), z) + 3/[o,oo)(^)- (3) 

The energy dissipation density due to damage progression is given by p where we assume the 
structure 


p(i) = -az + ^|zp -f /(-oo,o](i) 

with a,/3 > 0. Because of the quadratic term f |zp the damage evolution is called rate-dependent 
whereas /3 = 0 would correspond to rate-independent systems. See [Mie051lEM061IMR.Z101IMT10j 
for analytical results on rate-independent damage and numerical experiments (without phase 
separation). We also refer to [BSS051 IFK09| for further analytical investigations of damage 
models. In comparison to [HKllj we use a gradient-of-damage theory with the Laplacian Az in 
® instead of a p-Laplacian div(|Vz|^ ^Vz) with p > n. 

In conclusion, the systems we would like to consider in this work are governed by o, @ and 
m and can be rewritten as 

dtc = —Sw 

w = P(-div(rVc) -f Wfic) + IF®'(e(w), c, z)) 
div((T) = 0 

dp{dtz) - Az + W^(e(u), z) -h 5/[o,oo)(^:) 9 0 

where w denotes the chemical potential. Here, the matrix P denotes the orthogonal pro¬ 
jection of onto the tangent space TE = {x £ \ = O} of the affine plane 


in Hr,' 
in Dt’ 
in rix 
in Ht, , 


(^o) 
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S := {a: G | = l}- The operator S determines whether we have an Allen-Cahn 

or a Cahn-Hilliard type diffusion of the system. More precisely, 

Allen-Cahn: 5 : ^ R^), 5(/) := M/, 

(4) 

Cahn-Hilliard: 5 : R^))*, S{f) := (MV/, V -)^2 ■ 

In the Cahn-Hilliard case, the operator S is invertible when restricted to S : Y ^ 1) where the 
spaces Y and V are defined as 

{cGifi(H;R'^)| / c = 0, ^Cfc = 0}, 

k=l 

|c* G (iJ^(H;R^))* I {c*= 0 for all c = d(a;)(l,..., 1), 

where d G and for all c = Cfc, fc = 1,... 7v| Cfc : fc-th unit function. (5) 

We need to impose some restrictions on the mobility matrix M. We assume that M is symmetric 
and positive definite on the tangent space TE. In addition, due to the constraint = 1, 

M has to satisfy the property ^ki = 0 for all fc = 1,..., iV. Note, that M = MP. 

We abbreviate Dt '■= (0,T) x D and {dfl)T ■= (0,T) x dfl, where D C Oil with TdY‘~^{D) > 0 
denotes the Dirichlet boundary. The initial-boundary conditions of our systems are summarized 
as follows: 

c(0) = c° in H, 
z(0) = in H, 
u = b on Dt, 

and additionally for Cahn-Hilliard systems 

MVw • l/ = 0 on (3 H)t, (IBC) 

where 7^ is the unit normal on dfl pointing outward and b the boundary value function on the 
Dirichlet boundary D. The initial values are subject to 0 < zq < 1 and c° G E fl R^o a.e. in H. 
In the following, we assume that b can be suitably extended to a function on tdx- 

The paper is organized as follows: In Chapter 2, we introduce an appropriate notion of weak 
solutions for the system ( |5'oD . To handle the differential inclusion rigorously, we adapt the 
concept of energetic solutions originally introduced in the context of rate-independent systems 
(see for instance IMieP.^j l to phase separation systems coupled with rate-dependent damage. This 
approach was firstly presented in |HKII) . The main result and their assumptions are stated at 
the end of Chapter 2. 

In Chapter 3, we prove existence of weak solutions for a regularization of system 
in classical formulation as 

dtc = —Sw in 

w = P(-div(rVc) + Wf’P°\c) + W^^{e{u), c, z) + edtc) in 

div(cr) -I- edivdVupVu) = 0 in 

dp{dtz) -lYz- ediv{\Wz\P-'^Vz) + W;^^{e{u),z) + dl[o,oo){z) 3 0 in 


expressed 


Hr, 

Hr, 

Hr, 

^T, . 


(Se) 


a -1^ = 0 on (9H)r \ Dt, 
rVc-lt = 0 on (aH)T, (IBC) 

Vz -1^ = 0 on (9H)r 


Y := 

V := 
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where and hK®' satisfy certain polynomial growth conditions and p > n. The initial¬ 

boundary conditions are 

(lIBCIl with (cr -I- e|VupVit) • = 0 instead of ct • = 0. (IBCg) 

It turns out that the weak solutions of the regularized system have the following regularities: 
c € Vm G L^(f2T;and Vz G LP(r2T;R") (with p > n as above). 

They are constructed by adapting the approximation techniques developed in [HKll] . 

The limit problem e \ 0 for ( |5 'eD corresponding to ( |5'oD with is solved in Chapter 

4. The displacement field u obtained in this process has (fl; ]R")-regularity in the first instance. 
To establish existence results for chemical free energies of logarithmic type, we prove a higher 
integrability result for Vit in Chapter 5, which is based on some ideas of |EL911 IGarOOl IGarn5b) . 

Finally, Chapter 6 is devoted to logarithmic free energies for the concentration c. Following the 
approach in [GarOOl fCarOSb] . we use a suitable regularization with polynomial growth of 

the logarithmic free energy density to obtain a solution for ( |5'oD . Using this regularization, 

the chemical components Ck become strictly positive in the limit. 

The notation, we will use throughout this paper, is collected in the following. 


Spaces and sets. 

luy(u) 

Br{A) 

Qr{xo) 

{/ = 0 } 




standard Sobolev space 

functions of which are non-negative almost everywhere 

functions of which are non-positive almost everywhere 

functions of K") which vanish on D C in the sense of traces 

open neighbourhood of A C R" with thickness R 
open cube {a: G R" | ||x — xqHoo < R} 

zero set {a: G U | f{x) = 0 a.e.} of a function / G L^(f2) defined up to a 
set of measure 0 and defined uniquely if / G for p > n as 

lUi-P(O) co(n) 

the set (0, T) x 


Functions, operations and measures. 


[/] + 

Im 

Xm 

{9*J) 

dE 

p* 


non-negative part of /, i.e. max{0, /} 
indicator function of a subset M C X 
characteristic function of a subset M C X 

classical derivative of a function W with respect to the variable e 

dual pairing of g* G R”))* and / G ITi’’'(U;R”) 

generalized Clarke’s subdifferential of E 

Gateaux differential of E 

Sobolev critical exponent for n > p 


diam((5) diameter of a subset Q C R" 

Hausdorff measure of dimension n 
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Lebesgue measure of dimension n 


£" 


2 Existence theorem 


2.1 Weak formulation 


The weak notion, we will derive in this section for the doubly nonlinear differential inclusion 
occurring in ( |5'o| ), is inspired by the concept of energetic solutions for rate-independent systems 
(see for instance [Mien5] b In the rate-independent setting, the differential inclusion is formu¬ 
lated by a global stability condition and an energy inequality. In [HKllj . we have introduced 
an approach which uses an energy inequality and a variational inequality to handle the rate- 
dependence coming from the viscosity term in the damage dissipation density function 

P- 

The corresponding Gateaux-differentiable free energy £ : x x fi 

—?► R and dissipation functional TZ : —>• R to system ( [5^ are given by 

£{u,c,z) := f irVc : Vc-I--I-IF'^*'(c)-f IT®\e, c, z) dx, 

In 2 2 

Ti{z) ■■= -az + ^|i|^ dx, 

with viscosity constants a,/3>0. To account for the constraints of z, we extend the functionals 
£ and 7Z above for analytical reasons by indicator functions: 


£{u, c, z) := £{u, c, z) + 


'[0,c 


3)(z) dx. 


7^(i) := 7^(i)-f / /(_oo,o] (^) dx. 


n 


n 


If we equip the space H^{£l)r\L^{£l) with the norm the generalized 

subdifferential d^^£ at a point {u,c,z) G il^(r2;R") x x (H^(n) nL°^(fl)) is 


af £■(«, c, z) = {dj{u, c,z)+r£ {H^n) n 


€ dl. 


i/i(n)nL”(n 


)W}- 


( 6 ) 


The inclusion L^{£1) C fl L°°{H))* will be later used for the construction of a specific 

subgradient. Using property the differential inclusion in ( |5'oD can be rewritten in a weaker 
form as 


0 G d^^£{u{t), c(t),z(t)) + diTZ(z{t)). 

The analytical basis for a formulation of a weak solution is the following proposition (a proof of 
a related result can be found in |HKllj h 

Proposition 2.1 Let (u,c,w,z) G C^(Ot;R"' x R^ x R^ x R) be a smooth solution satisfying 
0 and @ with the initial-boundary conditions (lIBCI) . Then the following two conditions are 
equivalent: 

(i) 0 G d^^£{u{t), c(t), z{t)) dz'R-{z{t)) for all t G [0,T], 
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(ii) the energy inequality 


£{u{t),c{t),z{t)) + / {di'R{dtz),dtz)As+ / (5w(s),w(s))ds 
Jo Jo 

< f ('u(O), c(0), z(0)) + / W^^{e{u),c, z) : e{dtb) dxds 

J Ot 

for all 0 < t < T and the variational inequality 

0 < (dz£{u{t),c{t), z{t)) + r(t) + di7^(atz(i)), 
for all C € H}_{n) fl L°°(n) and r(t) € dlHf{Q)nL°°{n)izit)) for all 0 < t < T. 
If one of the two conditions holds then the following energy balance equation is satisfied: 
£{u{t),c{t),z{t)) + [ {diii{dtz),dtz) ds + [ {Sw{s),w{s))ds 


= £{u{0),c{0),z{0)) 


W^^(e(u),c, z) : e(dtb) dxds. 


Remarks for Proposition [2711 In contrast to |HK11] . the energy inequality in (ii) compares 
the energy at the beginning s = 0 with the energy at an arbitrary time s = t instead of s = 
with s = ^2 for 0 < < ^2 < T. 

Applying the chain rule on the right hand side of 


£iu{t), c{t),z{t)) - £(m(0), c(0), 2 :( 0 )) = 


dt 


£{u{s), c(s), z(s)) ds 


and using m, m as well as the variational inequality in (ii), the “>”-part of the energy balance 
can be shown. 


We will see that in our approach the mathematical analysis of (j^J) requires several e-regularization 
terms (see ( |A'eD ) to establish the energy and variational inequality for the differential inclusion 
and to handle the logarithmic free energy. A transition to e \ 0 will finally give us a solution of 
the limit problem diSpI). 


Proposition l2.1l can also be formulated for the regularized system (j^j) with the regularized energy 

£e{u,c,z) := [ ipvc : Vc+l\Vzf + W'=‘'’P°^(c) -b W‘'\e,c,z) + ||Vw|^ -b -jV^l^dx, 

Jo 2 2 4 p 

£e{u,c,z) := £e{u,c,z) + / /[0,oo) (2:) dx, 

Jn 

and the initial-boundary conditions dlBCgl ). Notice that Fdtc = dtc because of dtc{t,x) G TE. 

We can now give a weak notion of ( [671 ) and ([^j). (The energy densities and W®' will 

satisfy some polynomial growth conditions which are specified in the next subsection.) 

Definition 2.2 (Weak solution for the regularized system d6'eD ) We call a quadruple q = 
{u,c,w,z) a weak solution of the regularized system d6eD with the initial-boundary conditions 
dIBCeD if the following properties are satisfied: 
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(i) the components of q are in the following spaces: 


u&L°°{Q,T-W^'\n-W^)), u\D^=b\D^, 

ce L~(0,T;77i(0;]R^))ni/\0,T;L2(f];R^)), c(0) = c°, cGS a.e. in 
z e L^{0,T;Wl’^{n)) n H\0,T;L^{n)), z{0)=z°, dtz<0, 


and 


w G 


for C-H systems, 
for A-C systems 


(ii) for all C, G R'^) and for a.e. t G [0, T]: 


dtcif) ■ C, da; 


J^M'Vw{t) : VC da; 
Mw(t) ■ C dx 


for C-H systems, 
for A-H systems 


(Hi) for all f G and for a.e. t G [0,T]; 


w{t) ■ C da; 


PrVc(t) : VC + PVF,'i^’P°‘(c(i)) • Cda; 


+ [ ¥W)){e{u{t)), c{t), z{t)) ■ C + edtc{t) ■ ( da; 
JQ 


(7) 


( 8 ) 


(iv) for all ( G for a.e. t G [0,r]; 


/ W'^){e{u{t)),c{t), z{t)) : e(C) + e\Wu{t)\'^Vu{t) : VC da; = 0 

Jq 

(v) for all C G W^’^(r2) and for a.e. t G [0,T]; 

[ {e\\7z{t)\P~'^ + l)V2(t) • VC + {W^){e{u{t)),c{t), z{t)) - a + f3{dtz{t)))( dx 
JQ 


(9) 


>-{r{t),0, (10) 

where r{t) G (VF^’P(r2))* satisfies {r(t),z(t) — C) >0 for all C G W')J:’^(r2) 

(vi) energy inequality for a.e. t G [0,r]; 

£e{u{t),c{t), z{f)) - £e{u°, c°, z°) + [ a{z° - z{t)) dx 

Jn 

+ f /3|i9t2p + e|i9tcp da;ds + f {Sw{s),w{s)) ds 
Jilt Jo 

— [ W)){e{u),c, z) : e{dtb) dxdss ( |VupVM : Vi9t6da;ds, (11) 

J J 

where is the unique minimizer of £e{-,c^,z^) in with trace u^\d = 6(0)|n- 


With the help of the operator S, the diffusion equation © can also be written as 

[ dtc(t) ■ Cdx = 

JQ 

which will be used in the following. 

Definition 2.3 (Weak solution for the limit system ( |5'oD ) A quadruple q = {u,c,w,z) is 
called a weak solution of the system ( |5*oD with the initial-boundary eonditions (lIBCIl if the fol¬ 
lowing properties are satisfied: 

(i) the components of q are in the following spaces: 

cG L°°(0,T;iJ^(n;R")), cG E a.e. m fir, 
z e L°°{0,T;Hl{n))r\H\0,T;L\n)), z{0)=z°, dtz < 0 

and 

w G for C-H systems, 

w G for A-C systems 

(ii) for all C G with dtC G and ({T) = 0; 

f {c — c^) ■ dtC dx dt = f {SwX)dt 

J Ot’ J 0 

(Hi) for all ( G fl and for a.e. t G [0,T]; 

[ wit) • C dx = [ PrVc(t) : VC + VWficit)) ■ f dx 
Jn Ja 

+ [ FW(){e{u{t)),cit),zit)) ■ (dx 
Jn 

(iv) for all ( G and for a.e. t G [0,T]; 

f W()ieiu{t)),c{t),z{t)) : e(C) dx = 0 
Jn 

(v) for all ( G Hf{n) fl L°°{n) and for a.e. t G [0,T]; 

[ Vz(t) • VC + (W®J(e(u(t)), c{t), zit)) - a + P{dtzit)))( dx > -{r{t), (), 

Jn 

where r(t) G (i/^(0) fl L°°(0))* satisfies (x(t), z(t) — C) > 0 for all ( G fl L°°(n) 

(vi) energy inequality for a.e. t G [0,r]; 

£{u{t),cit), z{t))-\- f a{z^ — zft)) dx-\- f I3\dtz\'^ dxdsf {Swis),w{s)) ds 
Jn Jnt Jo 

< f (m°, c°, 2 ;°) + f W))ie{u),c,z) ■. e{dtb)dxds, 

J Qj 

where is the unique minimizer of £[■ ,c^, z'^) in H^i£l-,W^) with trace = 6(0 ) |d. 


9 




Note that both notions of weak solution imply chemical mass conservation, i.e. 


/ c{t) dx = const. 

Jn 

2.2 Assumptions and main results 

The general setting, the growth assumptions and the assumptions on the coefficient tensors which 
are mandatory for the existence theorems are summarized below. 

(i) Setting 


Space dimension 
Components in the alloy 
Regularization exponent 
Viscosity factors 

Domain 

n G N, 

iV G N with TV > 2, 

p> n, 

a,f3 > 0, 

D C M" bounded Lipschitz domain. 

Dirichlet boundary 

Time interval 

DC an with > o, 

[0, T] with T > 0 

(ii) Energy densities 


Elastic energy density 

W®' G X X R;M+) with 

W<^\e,c,z) = W<^\e\c,z), (Al) 

lV'='(e,c,z)<C(|ep + |c|2 + l), (A2) 

??|ei - e2p < (lTe(ei,c,z) - VE_®J(e2, c, z)) : (d - 62), (A3) 
|lT_^i(ei+e2,c,z)| < c, z) + |e2| + 1), (A4) 

|lT,^i(e,c,z)|<C(|ep + |cp + l), (A5) 

|lT,^i(e,c,z)|<C(|e|2 + |c|2 + l), (A6) 


Chemical energy densities G C\R^]R) with > -C, 



\Wf'^°\c)\<C{\cfl^ + l), (A7) 

lych.log(^) = 0 ^ Cfe log Cfc + -c • Ac, d > 0, A G K^ym” (AS) 

k^l 

(iii) Tensors 


Mobility tensor 

M G symmetric and positive definite on TS and 

Energy gradient tensor 

N 

Mki = 0 for all = 1 ,..., A^, 

r G £(R'^^"; R'^^") symmetric and positive definite 
fourth order tensor 
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Remark 2.4 Due to the effeet of damage on the elastic response of the material, Vk®* is often 
modelled by the following ansatz: 

where : [0,1] M+ is a continuously differentiable and monotonically increasing function with 
$( 0 ) = 0 and fj > 0 is a small value. 

A typically form of the elastically stored energy density Vk®* is as follows: 

W^\c,e) = ^{e-e*{c)):C{c){e-e*{c)). (12) 

Here, e*(c) denotes the eigenstrain, which is usually linear in c, and C(c) € £(K."y^) is a fourth 
order stiffness tensor, which is symmetric and positive definite. The elastic energy density is 
called homogeneous if the stiffness tensor does not depend on the concentration, i. e. C(c) = C. 

Note that the inhomogeneous elastic energy dni) fits into our setting with the previous growth as¬ 
sumptions (HB-dMl)- In particular, we are not confined to homogeneous elasticity as in \HK1 If . 
There, the more restrictive growth condition \ W^f{e,c,z)\ < C{\e\ + |cp + 1) is used instead of 

(KSI) . 


The main results of this work are summarized in the following theorems: 

Theorem 2.5 (Existence theorem - polynomial case) Let the above assumptions be satis¬ 
fied. Then for every 

c° G H^(f2;R'^) with c° € E a.e. in LI, 

£ H^{L1) with 0 < z^ < 1 a.e. in LI, 

there exists a weak solution q of the system ( | 6 'oD with and the initial-boundary 

conditions (lIBCI) in the sense of Definition \2.3[ 

Theorem 2.6 (Existence theorem - logarithmic case) Let the above assumptions be satis¬ 
fied and, additionally, let D = dLl and T = 7 ld with a constant 7 > 0. Then for every 

bGW^'^{Q,T-,W^'°°{Ll-W)), 

(P £ H^{Ll-, K'^) with c° £ E and c° > 0 a.e. in O for k = 1,..., N, 

z^ £ H^{Lt) with 0 < z^ < 1 a.e. in LI, 

there exists a weak solution q of the system diSpI ) with VF®*' = and the initial-boundary 

conditions jmsi) in the sense of Definition \2.tA Additionally, Cfe > 0 a.e. in LIt for k = 1,..., N. 

Remark 2.7 Note that for Theorem \2.6\ the assumptions (IA2I) . (IA5I) and (IA6II can be replaced 
by 

lT®'(e,c,z)<C(|ep + l), (A2’) 

|lT®'(e,c,z)|<C(|ep + l), (A5’) 

\Wfl{e,c,z)\<C{\e\^ + l), (A6’) 

for all c £ with 0 < Cfe < 1 and ^ ^ all z £ R with 0 < z < 1. 
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3 Existence of weak solutions of ([^ 


The proof is based on [HKllj . Arguments similar to |HK11) are only sketched. 

Since e > 0 is fixed in this section, we omit the e-dependence in the notation, e.g. £ always 
means here £e and so on. Furthermore, is assumed to be in W^'P{£1) in this section. 


1. Step: constructing time-discrete solutions. 

Set to be a minimizer of u £(u,c°,z^) defined on the space W^''^{VL) with the con¬ 
straint u\d = ^(0)|d in the sense of traces. 

Let the closed subspace of x x be defined by: 


u € 

QZ = { 

z e VFi’P(n) 


u\d = b[mT)\D, 

J^c — dx = 0 for C-H systems. 


0 < 2 < 


Based on the initial triple (u°, c°, 2 °), we construct (u^, d^,z'^) for to = 1,..., M recur- 
sively by minimizing the following functional —>■ M: 


E^(?x, c, z) := 8{u,c, z) + rH 


z — z 


M 


c — c 


M 


X 


£T 

T 


c — c 


M 


(13) 


L2 


where X denotes the space V (see ®) with the scalar-product 

(ci I C 2 )js: := / MV5“^ci • V5“^C2 da; 

Jn 

for Cahn-Hilliard systems and X = L^(n;R^) with the scalar-product 


(ci \c2)x ■= / Mci • C 2 da; 

Jn 

for Allen-Cahn systems. 

Note that the last regularization term in (1131) is not necessary for Allen-Cahn equations 
due to the term with the X-norm. To use a uniform approach, we consider this term in 
both systems. By direct methods of calculus of variations the triple 


{um,Cm,zZ)-= argmin E^{u,c,z) 

(u,c,z)eQ^ 


exists, cf. IHKIIj . Furthermore, we set 


Wm ■= 


-S 


( m m — 1 \ 


+ A^ 


M 


-1 


for C-H systems, 
for A-C systems. 


with the Lagrange multiplier (associated with the mass constraint for C-H systems) 
given by 

A)& := / + W:^Mn^),CM,zZ)dx. 

Jn 
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We define the time incremental solutions as 

and introduce the piecewise constant interpolations qm, Qm-, tM, t~M linear interpo¬ 

lation qM as 

Im '■= min{TOr | m G No and mr > t}, 

:= min{(m — 1)t | m € Nq and mr > t}, 
qiviit) := qZ for t € {{m - l)r,mT], 

9m(^) := 9m for t & [^t, (m -P 1 )t), 

9M(t) := PqZ + (1 - /3)9m“^ for tG [(m - l)r,mr) and /3 = ^ - (to - 1). 


Due to the minimization properties of {u'^, z'^), we establish the following variational 

formulas and energy estimate (cf. |HK11[ Lemma 6.2]): 


Lemma 3.1 (Euler-Lagrange equation, energy estimate) The functions qM, qM “rid 
qM satisfy the following properties for all t G (0,T): 


(14) 


fij for all ( G 

[ {dtCM{t)) • Cdx = -{SwM{t),0 
JQ 

(li) for all C e 

[ WM{t)-Cdx= [ PrVcM(t) : VC + PtE_"^P°HcM(t))-Cdx 

Jo, Jq 

+ / PW®c(e(MM(t)),CM(t),^M(t)) • C + e5tCM(t) ■ Cdx (15) 
Jn 

(Hi) for all f G 

[ W^){e{uM(t)),CM(t), ZM{t)) : e{C) + e\VuM{t)\‘^'^UM(t) :VCdx = 0 (16) 

Jn 

(iv) for all ( G W^’P(fl) with 0 < C -I- ZM{t) < z)(j{t): 

[ {e\VZMit)\^~‘^ + l)VzM(t) • VC + W^){e{uMit)), CM(t), ZM(t))C dx 

Jn 

+ [ {-a + l3{dtZMit)))Cdx >0 (17) 

Jn 

(v) energy estimate: 

r*M r (3 £ 

£{uM{t),CM{t),ZM{t)) + / -adtZM + ^\dtZM\^ +-zldtCMl^ dxds 

Jo Jn z z 

ptM ^ 

+ J -{SwM{s):WM{s))ds-£{u°,c°,z°) 
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< 


piM n 

/ / ^^,e'(e(MM + ^ - ^m)> Cm: ^m) : e{dtb) dxds 

Jo Jr2 

ptM r 

+ e / |Vu^ + V6-V&^pV(u;;;;f + 6-6;;^) : V9t&dxds. (18) 
Jo Ja 


2. Step: identifying convergent subsequences. 

The energy estimate (v) in Lemma 13.11 growth condition (IA4I) and a Gronwall estima¬ 
tion argument lead to a-priori estimates for the energy £{uM{t),CM{t), ZM^t)) and for 
||9tZM||L2(ni,), WdtCMWh'^iQ.T) and {Swm{s),wm{s)) ds. By sta ndard compactness argu¬ 
ments and a compactness theorem from Aubin and Lions [Sim86) , we deduce the following 
weak convergence properties, cf. |HK11) : 


Lemma 3.2 There exists a subsequence {M^} and an element q = {u^c,w,z) satisfying 
(i) from Definition \ 2 .‘A such that for a.e. t G [ 0 ,T]; 


(t) UM,^u in L°°{ 0 ,T;W^J{n)), 
^ c{t) in 

CMk ^ c in iL^( 0 , T; R^)), 


(zti) ZMk,ZM^ ^ z zn L^{0,T-W^’P{n)), 
ZMk(t),ZM^(t) zit) in W^'Pin), 
^Mk,ZM^ ->• z a.e. in ftr, 

ZMk ^z in H\ 0 ,T;L\n)) 


(iv) iCMfc ^ IT m L^( 0 , T; R^)) 

WMk ^ w in L^(fl 7 ’;R^) 
as fc —>■ oo. 


for C-H systems, 
for A-C systems 


Exploiting the Euler-Lagrange equations, we can even prove stronger convergence proper¬ 
ties. To proceed, we recall an approximation lemma from |HK11) . 

Lemma 3.3 f |HKlll Lemma 5 . 2 ]) Let q > 1 , p > n and /, C S L'?( 0 , T; with 

{C = 0 } 2 {/ = 0 }. Furthermore, let {/m}mgn C L'J( 0 , T; LEj;’^(f 2 )) be a sequence with 
fAiif) f(t) in W^’P(fX) as M ^ oo for a.e. t € [ 0 ,T]. Then there exists a sequence 
{Cm}mgn C L'^( 0 , T; LEj:’^(f 2 )) and constants VM,t > 0 such that 

(i) Cm in L‘>{ 0 , T; W^’P{n)) as M ^ oo, 

(a) Cm < C n.e. in fix for all M gN, 

(Hi) VM,tCM{t) < fnif) a.e. in fl for a.e. t G [ 0 ,T] and for all M gN. 

If, in addition, C f a.e. in LIt then condition (Hi) can be refined to 
(Hi)’ Cm < fM a.e. in LIt for all M gN. 

We are now able to prove strong convergence results by using uniform convexity estimates. 


Lemma 3.4 (Strong convergence of the time incremental solutions) There exists 
a subsequence {Mk\ such that for a.e. t G [ 0 ,T]; 
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UMk{i),UM^{t) u{t) in 

“-s- ^T, 

(a) CMk,CM^ -J> c m L2*(0,T;i?i(rj;R^)), 
CMfe(t),c^Ji)c(i) in 
CMfc,c^j, —>■ c a.e. m fir, 

CMfc ^ c in H^{0,T; L‘^{Vl\ R^)), 
as k ^ oo. 


(m) Zm„Zm^ ^ z m LPiO,T;W^’Pin)), 

w W^'P{Vt), 
z a.e. in Qt, 

ZM, ^z in H\0,T-L^{n)) 


Proof. We omit the index k in the proof. 

(i) We refer to |HK111 Lemma 5.9]. 

(ii) The weak convergence properties for cm, and Cm^ follow from Lemma [3.21 It 
remains to show strong convergence of Vcm to Vc in L^(r 27 ';R^). 

By the compact embedding H^(D,]'R^) ^ /^+^(r2;R^) and Lemma [3.21 we get 

||CM(t) - c(t)|| j^ 2 */ 2 +i(f 2 .][ 5 N) —>■ 0 as M —>• oo for a.e. t € [0,T]. The boundedness 
property esssup^gjQ .j.]||cM(t) — c(t)||i 2 */ 2 +i(Q.RAr) < C for all M G N and Lebesgue’s 
convergence theorem yield cm —t c as M —>• oo in Testing (ITS]) with 

C = Cm (t) and with ^ = c{t) gives after integration from t = 0 to t = T: 


/ PPVcm : Vcm dxdt = / wm ■ cm — PIT^*'’P°^(cm) • cm dxdt 

Jut JtlT 

- / PW^(e(uM), cm,zm) ■ CM + edtCM ■ cm dxdt, 
J Ot-i 

j PPVcm : Vcdxdt = f wm • c — PW^^’P°*(cm) • cdxdt 

J J Cl-j’ 

- / PW^(e(uM), cm,zm) ■ C + edtCM ■ cdxdt. 

J O7-' 

Passing to M —>■ 00 and comparing the right sides of the equations shows 


PPVcm : Vcm dxdt — t / PPVc : Vcdxdt. 


By using the properties PVcm = Vcm and PVc = Vc, we eventually obtain 


rVcM : Vcm dx dt —t 


rVc : Vcdxdt. 


We end up with 



r(VcM — Vc) : (Vcm — Vc) da: dt —>■ 0. 


Therefore Vcm —t Vc in L^(fl'r;R^) since P is positive definite. 

(iii) Applying Lemma 13.31 with f = z and /m = z]^ and C = z gives an approximation 
sequence {Cm} ^ L^’(0, T; W:J:’^(fI)) with the properties: 


CM^^in LPiO,T;W^’P{n)), 


(19a) 
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0 < Cm < for all M €N. 


(19b) 


The estimate 

Cucl^ZM - Vzr < (|VzMr-"VzM - |Vz|P-2Vz) • V(ZM - z) 

where Cue > 0 is a constant and equation (HZl) tested with C = Cm (t) — zm (t) (possible 
due to (I19bl) l yield: 

Cue / ejV^M — dxdt + / dxdt 

< f ((eIV^mT”^ + 1)V2m — (£|V2;|^“^ + l)Vz) • V(zm — 2:) dxdt 

< f (elV^M^”^ + l)VzM ■ V(2 :m - Cm) dxdt 

J Ot' 

+ [ {e\\7zM\^~^ + 1)VZM • V(Cm - z) - {e\Vz\P-^ + l)Vz • V(zm - z) dxdt 

J 0^1 

< / (W^®i(e(MM),CM,^M) - a +;d9t£M)(CM - 2M)dxdt 

+ [ {e\\7zM\^~^ + 1)VZM • V(Cm - z) - (£|Vz|p- 2 + l)Vz • V(zm - ^) dxdt 

J 0^1 

< ||lT®^(e(MM),C m, Zm) — a + l3dtZM\\L'^(nT) IICm — 2m||l2(q^) 


bounded 


+ (e||VzM||^p(Qy) + \\^ZM\\Lv/iv-i)(nT)) II^Cm - V2 :||lp(Oj 


bounded 


f (e|Vz|^ ^ + l)Vz • V(zm —- 2 ) dxdt 

J nj" 


Due to (I19al) and zm ^ z in L°°(0, T; bb^’P(D)) as well as zm ^ z in L^(Dr), each 
term on the right hand side converges to 0 as M —>■ oo. ■ 


3. Step: establishing a precise energy inequality. 

In this step we establish an asymptotic energy inequality, which is sharper than the energy 
inequality in (fT^ . Note, that compared to (fT^ the factor 1/2 in front of {Swm{s), wm{s)) 
is missing. To simplify notation, we omit the index k in the following. 


Lemma 3.5 For every t € [0,r]; 

ptM 


pZM P 

£{ uM { t ), CM { t ), ZM { t )) + / / -adtZM + PldtZM ^ + eldtCM^ dxds 

Jo Jn 

ptM 

+ {SwMis),WM{s)) ds - £{u°,c°,z°) 

Jo 
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< 


ptM r 

/ / + ^ - ^m)> Cm: ^m) : e{dtb) dxds 

Jo Jq 

ptM r 

+ s / |Vm^ + V6 - V6^pV(u^ + 6 - 6^) : V^t^dxds + Km 
do Jq 


lo JQ 
with Km —>■ 0 as M —>■ oo. 

Proof. Apply ing the estimate lE^(gM) ^ “ ^M~^CM>^M) for m = 1 to 

^ yields (cf. [HKlll Lemma 6.10]): 

,0 „o .,0\ 


£{uM{t),CM{t),ZM{t)) - £{u^,c^,z^) 
ptM 

< s 

h Jfi 


pT^M P 

/ / NiuM + bis)-blj)\'^V{u^ + b{s)-blj) :Vdtb{s)dxds 

Jo Jn 

ptM P 

+ / ^.e (e(MM + b- ^m): Cm, ^m) : 6(9* *6) dxds 

70 -/Q 

ptM P 

+ / / lL,f(e(a)(, + 6 m - 6j^^), CM, -Zm) • ^tCM dxds 

do dn 


(*)i 


ptM P 

/ / rVcM : ydtCM + W",c^’‘’°'(cm) • dtCM dxds 

Jo Jo. 


(*)2 


ptM 


'O JO 


W^l{e{Ufj + bM - bj^j),CM, ZM)dtZM dxds 


(**)i 


+ 


ptM P 

/ / ffjVzMr^^VzM • VStZM + Vzm ■ V9tZM dxds . 

do Jn 


( 20 ) 


(**)2 


The elementary inequalities 

(|ViMr“^VfM - IV^mT'^V^m) ■ VStZM < 0 and (Vzm - ^zm) dtZu < d 
and (fT71) tested with ^ := —dtZM{t)T lead to the estimate: 

(**)i + (**)2 


P^M P 

- / -adtZM + PldiZMl"^ dxds 

do Jn 

ptM P 

+ / + 6m - b'lj),CM, zm) - (e(MM), cm, ZM))dtZM dxds . 

Jo Jo ' 


Furthermore 

ptM 


pZM P 

(*)i < / / W^l{e{uM),CM,ZM)-dtCMdxds 

Jo Jo 

ptM P 

+ / {W‘'^{e{u~^ + bM-b]^),CM,ZM)-W^^{e{uM),CM,ZM))-dtCMdxds. 

Jo Jn ^ 
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Using the elementary estimate r(VcM — Vcm) : VdtCM < 0 gives 

rtM 


r^M r 

W2 < / / rVcM : V^tCM + ■ dtCM dxds 

Jo Jn 

rtM f 

+ / / {Wf’P°\cM) - Wf’P°\cM)) ■ dtCM dxds . 

Jo Jn ' 


Hence, applying equations (ITSl) with ( = dtCM{t) and (|n)) with ( = WM^t) by noticing 
fdtCM{t) = dtCM{t) shows 

rtM rtM r 

{*)i + {*)2<- {SwM{s),WM{s))ds- / e|5tCMp dxds + + Km- 

Jo Jo Jn 

Lebesgue’s generalized convergence theorem, growth conditions (IA5ll - (IA7l) and Lemma[3Tl] 
show Km := + ^^m + 0 as M —>• 00. We would like to emphasize that we need the 

boundedness of Vmm in and the boundedness of dtCM and dtZM in 

with respect to M. ■ 


4. Step: passing to M — >■ 00. Using Lemma [3.21 Lemma [3.41 and dT41) . (TTSl) and (ITSl) we 
establish (ii), (hi) and (iv) of Definition o Moreover, Lemma 13.51 implies 


£{uM{t),CM{t),ZM{t)) + / -adtZM + fJ\dtZM\'^ + £\dtCM\'^ dxds 

J 0,-1- 

+ [ {SwM{s)jWM{s))ds-£{u^,c^,z^) 

Jo 

rtM r 

< / / (e('UM + b- 6m), Cm,zm) ■ e{dtb) dxds 

Jo Jo 

rtM r 

+ e / |Vmm + V6-V6mPV(um + 6-6m) : V3t6dxds + Km- 

Jo Jn 


The energy estimate (vi) from Definition 12.21 follows from above by using the known con¬ 
vergence properties and weakly semi-continuity arguments. 

It remains to show (v) of Definition 12.21 To proceed, we cite the following lemma from 
[HKll] which provides a tool to drop a restriction on the space of test-functions for a 
variational inequality of a specific form. 


Lemma 3.6 ( |HKlll Lemma 5.3]) Let p > n and f € Lp^^p ^)(D;R"'), g € L^{Q), 
z € lUb^(D) with z > 0, f ■ Vz > 0 and {/ = 0} D {z = 0} a.e.. Furthermore, we assume 
that 

[ f ■ + 3 C da: > 0 for all ( € lUl’^(D) with {(^ = 0} D {z = 0}. 

Jn 


Then 


//•VC + 5 Cdx> / bj+Cda: for all f e (fl). 

Jo J { 2 —0} 
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We are now able to prove the remaining property. 


Lemma 3.7 We have 

[ {e\\7z{t)\P~'^ + l)Vz{t) ■ VC + (VL®i(e('u(t)), c(t), z{t)) -a + f3{dtz{t)))C dx 
JQ 

>-{r{t),0, (21) 

for all C £ and for a.e. t G [0,r], where r{t) £ L^{Q) C is given by 

rit) := -X{z{t)=o}[Wf^{e{u(t)),c(t),z(t))] + . (22) 


Proof. First of all, we take any test-function C £ LP{0,T;Wb^{fl)) with 

{C = 0} 2 {-2 = 0}. Lemma [S31 gives a sequence {Cm} C LP{0,T;W^’^{n)) with Cm ^ C 
in LP{0,T;W^'P{il)) and 0 > vCmit) > —ZM{t) where v depends on M and t. Therefore 
(HZD holds for C = CM{t)- Integration from 0 to T and passing to M —>• oo gives 

[ (e|V 2 |P -2 + i)Vz . VC + {W^Xe{u), c,z)-a + P{dtz))C dxdt > 0. 

J 

In other words, 

[ (e| V 2 (t) r-2 + l)Vz{t) ■ VC + Wf^{e{u{t)), c{t), z{t))C dx 
Jn 

+ f {-a + l3{dtz{t)))Cdx>0 
Jq 

holds for every C £ with (C = 0} 3 {z{t) = 0} and a.e. t £ [0,r]. To finish the 

proof, we need to extend the variational inequality to the whole space Wb^{VL). 

Setting / = {£\Vz{t)\P~'^ + l)Vz(t) and g = Wf^{e{u{t)), c(t), z(t)) — a + f3{dtz{t)), Lemma 
13.61 shows for every C £ IF}:’^(0) 


{e\yz{t)\P ^ -f l)Vz(t) • VC + {Wf^{e{u{t)),c{t),z{t)) - a + P{dtz{t)))C dx 

> [ [Wf^{e{u{t)),c{t),z{t)) - a + ^{dtz{t))]+Cdx 

i{z(t)=0} 

> [ [W^^{e{u{t)),c{t),z{t))] + Cdx. 

J{zit)=0} 


Now, variational inequality m follows by setting 

r(t) := -X{z{t)=o}[W^zieiu(t)),c{t), z{t))] +. 


Remark 3.8 Lemma \3. 7| gives more information than (v) from Definition \2.‘A It provides 
a special choice for r[f) given by 
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4 Existence of weak solutions of ([S^ - polynomial case 


In this chapter, we show that an appropriate subsequence of the regularized solutions for 
e G (0,1] of Definition 12.21 converges in “some sense” to q which satisfies the limit equations 
given in Definition [531 Besides that the initial damage profile in this chapter is in We 

approximate 2 :° G by a sequence { 2 °} in such that 2 ° —>■ 2 ° in H^{n) as e \ 0. 

Using the energy inequality and Gronwall’s inequality, we establish again the following energy 
estimate. 

Lemma 4.1 IVe have 

£e{Ue{t),Ce{t),Ze(t)) + f f -adtZe + l3\dtZe\‘^ + e\dtCe\'^ dx ds + f {SWeis),Weis))ds 

Jo Jq Jo 

<c(£:,(«°,c°,2°) + i) 

for a.e. t G [0,T] and every e G (0,1]. 

Since c°, 2 °) < £e{ui,c°, z^) < c°, 2 °), the left hand side is also uniformly bounded 

with respect to a.e. t G [0, T] and every e G (0,1]. By using standard compactness theorems and 
uniform convexity properties of lU®' (see issi)), we obtain the following convergence properties 
(cf. IHKllj h 

Lemma 4.2 (Convergence properties of q^) There exists a subsequence {sk} with Sk \ 0 
as k ^ 00 and an element q = (u,c,w,z) satisfying (i) of Definition such that for a.e. 

tG[0,T] 


(i) Uej, u in L^{0, T ; K”)), 

(iii) z,, ^2 inL°°{0,T-,H\f})), 

^Om L^{0,T-,L^in-,R^)), 

r-^V2e, ^Oin L^{0,T;LP{n)), 

Uef,{t) u{f) in H^{D;R^), 


Ugf. — >■ M a.e. in £1 t, 

Zef, —>-2 a.e. in fix, 

— >• in H^{Q;R^), 

^ 0 ot L‘^{n;R^), 

z,, ^zin H\0,T-L^{n)) 

(ii) J^cin L°“(0, T; ili(0;K^)), 



SkdtCet, 0 m 

Ce^if) cif) in 
Cef, c a.e. in £It, 


as k ^ 00. We additionally obtain for Cahn-Hilliard systems 

-^w in L‘^{ 0 ,T;H\n;R^)) 


and for Allen-Cahn systems 


Wek ^ w in L^iVlT^R^), 

^cin H\0,T;L‘^{n;R^)) 


as k ^ 00. 

As before, we will omit the index k in the subscripts below. 
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Remark 4.3 We would like to mention that the arguments in \HK11\ Lemma 6 .I 4 J cannot be 
adapted to prove strong convergence properties oJVce and due to the more generous growth 
condition (IA5I) as well as the use of Lemma \ 3. 3\ where the compact embedding ^ 

for p > n with a > 0 and a < 1 — j is exploited. 

We are now able to establish existence of weak solutions of ( | 5 'o| ) in the polynomial case. 

Proof of Theorem 12.51 Whenever we refer in the following to ©-(HID the functions u, c, w, z 
and r are substituted by Ue,Ci;,Wi:, z^ and r^. Moreover, Lemma 14.21 is used without mention in 
the following. 


(i) Let C S with dtf G and C{T) = 0. Integration from t = 0 

to t = T of © and integration by parts yield 


[ {cs - c°) ■ dtCdxds = [ {Sws,0ds. 
J Qt- J 0 


Passing to £ \ 0 shows (ii) of Definition 12.31 

(ii) Let C G fl Integration from t = 0tot = Tof© and 

passing to e \ 0 yield 


la 

Note that 


[ w(:dxds= [ PTVc: VC-H(PIT_"^‘'’P°Hc)+Pl^,c(e(M),c,z))-Cdxds. 

J Qt" QiT’ 


/ edtCe ■ C da;ds 

J ©T-- 


< £||9tC£||i2(Q.j,.RN)||^|H2(f2j,;RN) 0 

as £ \ 0. This shows (iii) of Definition 12.31 with 

(iii) Let f G ITj! 5 d(D;R") be arbitrary. Passing to £ \ 0 in ([HD yields for a.e. t G [0,T] 

[ W^ei4u{t)), c(t), z(t)) : e(C) dx = 0, 

Jn 

by noticing 


(23) 


£|VM£(t)pVue(t) : VCdcc 


< £||Vu£(<)|||4(f2)||CI|L4(n) 0. 


A density argument shows that (1231) also holds for all f G Therefore, (iv) of 

Definition 12.31 is shown. 

(iv) The characteristic functions X{z^=o} are bounded in L°°{LIt) with respect to £ G (0,1]. We 
select a subsequence such that X{ze = 0 } ^ X in i°°(flT) as fc —>■ oo. In the following, we 
will omit the index k in the notation. Integrating (HOD from t = 0 to t = T and passing to 
£ \ 0 show 



Vz-VC + (IT®He(u),c,z) 


a + f3{dtz))f dx > j x[W^.z*(e(M), c, z)]^C da:ds (24) 

J Cl-j’ 
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for all C e ^^(0, T ; l^l’^(fl)) fl L°°{U,t)- We also used the fact that 

[ e\yz,r^Vz,-yCdxds <£||Vz,||^;('^^)||VC|Up(o,)^ 0. 

J Qt-* 


It follows that 


Vz{t) ■ VC + (W®J(e(u(t)),c(t), 2 :(<)) - a + l3{dtz{t)))C,dx 

; 

/ x{t)[W^l{e{u{t)), c(t), z(f))] + C dx 
Jn 


> 


for all C S H^{n) n L°°{n) and a.e. t € [0,T]. Set r := —x[W^^{e{u),c,z)]^. For every 
C G F°°([0,T]) with C > 0 a.e. on [0,T] and every C S H^{n) fl L°°{D,) we also have 

0 > / ( [ - Zs{t))dx) ^{t)dt = f rs{(- ZE)^dxdt 

Jo \Jq, J «/Qt 

[ r{( — z)^dxdt = f f f r{t){( — z{t)) dx\ ^{t) dt. 

J Cl-j’ J 0 \J o J 

This shows r(t)(( — z(t)) dx < 0 for a.e. t G [0, T]. Hence, we obtain the inequalities (v) 
of Definition 12.dl 

(v) Weakly semi-continuity arguments lead to 

liminf ^£:£(ue(f),Ce(t), 2 :e(t))-I- f a\dtZe\ + l3\dtZe\'^ + e\dtCe\^ dxds + f {SWe,We)ds^ 

> S{u{t),c{t),z{t)) + f a\dtz\ + li\dtz\'^ + f {Sw,w)ds. 

iOt Jo 

Testing Q with (^ = — 6(0) and (iv) of Definition 12.31 with C = — 6(0) yield 

e f |Vit°|'‘da; = e [ |Vit°pVM° : V6(0) dx 
Jq Jq 

- [ w:,'(e(u°),c°,4 ):e(u°-6(0))dx 
Jn 

- f IF®‘(e(w°),c°,z°) :e(M°-6(0))dx = 0 

Jq 

as e\0. 

Therefore, we can pass to the limit e \ 0 in m and obtain (vi) from Definition nn ■ 


5 Higher integrability of the strain tensor 

To prove existence results for chemical free energies of logarithmic type, a higher integrability 
result for the strain tensor based on [GarOOl IGarOSb] will be established. We adapt the higher 
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integrability result for solutions of the elliptic equation of the form 

fdiv(W"_®'(e(M), c)) = 0 on fir, 1 

j Vb®i(e(M), c) • it = CT* • 1/ on {dn)T J 

to our setting with non-constant Dirichlet boundary data b and the additional damage variable 
z in ( |5'oD . In the following, we will use the assumption D = dfl. 

The proof of the higher integrability result is based on the following special cases of the Sobolev- 
Poincare inequalities and on a reverse Holder inequality. 

Theorem 5.1 (Sobolev-Poincare type inequalities) Let 1 < p < n. There exists a con¬ 
stant C > 0 such that 


(i) for all rectangles Q C R" and all u € W^’^{Q): 




(diamQ), 


(a) for all rectangles Q = rir=i(®»’^») — ^ ^ W^’P{Q) with u = 0 on 

{(xi,... , a„) \ at < Xi < bi, i = 1,..., n — 1} C dQ (in the sense of traces): 

(J^ WV ^ < C (diamQ). 

Theorem O can be obtained by considering the corresponding inequalities on the unit cube 
(0,1)" (for instance the case 1 < p < n was proven by Sobolev |Sob38] while Nirenberg [Nir59] 
gave a proof to p = 1) and then using a scaling argument. 


Theorem 5.2 (Reverse Holder inequality, see |Gia83| I LetQ C M" he a cube, g € 

for some q> 1 and g > 0- Suppose that there exist a constant b > 0 and a function f € Lioc(Q) 

with r > q and / > 0 such that 


/ 

JQr{xo) 


g'^ dx <h 



+ -f rdx 

'J Q2r{xq) 


for each xq € Q and all R > 0 with 2R < dist(xo, dQ). Then g € for s € [q,q e) with 

some e: > 0 and 


'Qr{xo) 





for all xo € Q and R > 0 such that Q 2 r{xq) C Q. The positive constants c,£>0 depend on b, 
q, n and r. 


Theorem 5.3 (Higher integrability) Let b € R"), z € L°°{Q) with 0 < z < 1 a.e. 

in H and c € L^{Tt;M.^) for some /i > 4. Then there exists some p € (2,p/2] such that for all 
u € IL^(Q;R") which satisfy u\d = b\n and 

[ IP"/(e(u),c,z):e(C)da; = 0/oranCeil/5(H;R"), (25) 

Jn 
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we obtain u G R”) and 

l|Vw||LP(0;R"Xn) < C(||VM||i2(Q.Rnxn) + 11 c| | ^2p(n;RW ) + !)■ (26) 

The positive constants p and C are independent of u, c, z. 

Proof. The proof is based on |Gar001 Lemma 4.4 and Theorem 4.3] and uses a covering argument. 
However, due to the non-constant boundary condition, we need to apply a more general Sobolev- 
Poincare inequality (see Theorem 15.II (iiP than in [GarOO] . 

(i) Higher integrability at the boundary. 

Let xq € dfl. Then there exist an i?o > 0 and a bi-Lipschitz function t : Q —)■ R” with the 
open cube Q := Qro(O) such that xq G t(Q) and 

r(Q+)CL!, 
r(Q-) CR-\n, 

where Q+ := {x G Q\xn > 0} and Q~ •.= {x G Q\xn < 0}. Define the transformed 
functions u,h G c G and z G L^{Q^) as 

(m,6,c,z)(x) := { u , b , c , z ){ T { x )). 

To proceed, let yo G Q and R < ^dist(2/o, and define for each R' > 0 the sets 

Q%iyo) ■■= {x G QR'iyo)\xn ^ 0}. 


We distinguish three cases: 

Case 1. We first consider the case (^^(j/o) 0 and Q~f^{yo) 0- 

The bi-Lipschitz continuity of r ensures 

dist(r(5Q+^(i/o)) n H, r(5Q^(?/o)) n D) > RCi, 

where Ci > 0 is independent of R and ?/o- Let f G C“(H) be a cutoff function with the 
properties: 

(a) ^ = 0 in H \ r(Q 2 R(yo)), (c) C = 1 in t(Qr(i/o)) G H, 

(b) 0 < C < 1 in H, (d) |V^| < ^R~^- 

Testing (OSl) with ( = — b), using the computation 

e(C) = ee(u) - ee{b) + ^((ii - b){WfY + - b)% 

and (lAll) . we obtain 


f (e(u), c, z) : e(u) dx 

Jq 

= [ f'^Wf^{e{u),c,z) : e{b)dx-2 f ^W®'(e(u), c, z) : ((m - 6)(VO*) dx. 
Jq Jq 

By iSSl), (EH) and (IMl) we also have the estimates 

v\e{u)Y < We(e(M),c,z) : e{u) + C{\cY + l)|e(M)|, 


(27) 
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IW^Mu), c, z) : iiu - 6)(ve)‘| < ^(|e(u)| + |c|^ + 1)|« - 6|, 
\W:^Mu),c,z):e{b)\<{\e{u)\ + \c\^ + lMb)\. 

Therefore, (EZl) can be estimated by 

V [ C^|e(w)pda; < C f + l)|e(u)| da; + ^ [ ^(|e(u)| + |cp + 1 )|m - 6| da; 

JQ Jci ^ Jn 

+ C [ C2(|e(«)| + |c|2 + l)|e(6)|dx. 
an 

Young’s inequality yields 

Cl [ dx < C [ + 1) da; + ^ [ \u- da;. (28) 

an an ^ an 

We choose /a = /q+ ^y^^udx. The calculation e{^{u — /a)) = ^e(ia) + ^((u — m)(V^)‘ + 
V^(ia — /a)‘) leads to 


/ |e(^(ia -/a))p dx < 2 ^ / f^|e(u)pdx+ / |u - /ap| dxV (29) 

./ o \J n Jq J 

Combining (1^ and (1^ . applying Korn’s inequality for i7^-functions with zero boundary 
values and using (a) and (b) gives 


/ |V(C(aa-Aa))pdx<C / 
an Jt( 


(|c|^ + l)dx+^ / ^ \u-b\^dx 


+ 


_c 

i?2 


•^(QJrIs/o)) 


; — /ap dx. 


Because of V(5(?a — /a)) = ^Via + (la — /a)(V^)* we derive by (a) and (c) the following type 
of Caccioppoli-inequality: 


>r{QUvo)) 


iVaapdx < C 


•^(Qiails/o)) 


+ — f 

aT(Q+j(yo)) 


(|c|^ + l)dx + ^ j^ 
\u — dx. 


. — dx 


Integral transformation by r implies 


[ |VMpdx<C'/' (I 

•'Qfiivo) •aQtofi/n) 


c 


Q^R^yo) 

C 


c,- + l)dx+^/ 




|zi — dx 


+ -^ I ^ \u-^^\^dx. 

R yQtR(vo) 


The condition Q 3 j^{yo) ^ 0 and D = i90 imply that u — b vanishes on i9(Qjj^(?/o)) CK" ^ x 
{0}. Therefore, we obtain by applying both variants of the Poincare-Sobolev inequality in 
Theorem 15.II for p = 2nl{n + 2): 

[ iVwpdx^C'/ (|c|^-hl)dx-h^/:’"(Qj^(yo))“"diam(Q+^(?/o))^ 

jQt(vo) JQtr,(vn) ^ 
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|Vu — V6| "+2 dx I + 




zn 

IVmI "+2 da; 


( 30 ) 


Note that if n = 1 we cannot apply Theorem 15.11 because of p = 2nl{n + 2) < 1. In this 
case, we can work with the inequalities in Theorem 15.11 where p is substituted by 1 and 
p* is substituted by 2. However, we will only treat the more delicate case n > 2 in the 
following. 

The estimates diam((52fi(2/o)) < CR and /1 ”(QJ/j(2/o)) > .R” (because of Q'^{yo) ^ 0) show 


''n(n+ A\^^(n+ 


'C”(QJr( 2 /o)) "diam(Q+j(yo)) <C. 
Now, dividing (1501) by C^i^Quiyo)) and using (1511) and 


(31) 


R2 >C"(Q2ij(2/o)) 


< C 


1 


^”(Q2fi(2/o)) 


gives 


^”(Qfl(yo)) JQ+^yo) 


[ iVfipda; < ^ ^ j (|c|^ + l)da; 

Jo±(vn) ^ \Q 2 Ryy 0 )) JQ+^iyo) 


Observe that 


C 


+ C 


1 


^C^iQ2Riyo)) jQ+^iyo) 
1 


C^iQ2Riyo)) JQ+j^ivo) 


I 

Jo: 


|Vm| "+2 da; 


I V&l "+2 da; 


\C^{Q2Riyo)) JQ+^iyo) 
Define the following functions on Q: 


|V5|^da; < 


g{x) := 


|Vit(a;)|"+2 for X € (3+, 

0 for X G Q \ (3+ 


and 




+2 for X G Q~^, 

for X G Q \ Q+. 


We eventually get 


/ dx < / 

Jqr(,vo) yQ 


dx + C 


QzRiyo) 


[i 


gdx 


QzRiVo) 


(32) 
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Case 2. Assume Q~^{yo) ^ 0 and Qsj^ivo) = 0- 
The bi-Lipschitz continuity of r implies 


dist(T(aQ|fl(yo)),T(aQ/i(yo))) > RCi, 

where Ci > 0 is independent of R and yo. Therefore, we can choose a cutoff function 
^ which satisfies 

(a) ^ = 0 in \ t(Q^j^(xo)), (c) ^ = 1 in T(Qji(xo)), 

(b) 0<C< 1 inf7, (d) |V^| < 

Testing (1251) with ^ — /i) and y ■= (xg) yields as in the previous case 

o ^ 


iVuPda; < C 


J t(Qr(xo)) 
Consequently, 


'^iQSfiixo)) 


(Id' 




r(Q3 R(a:o)) 


|m — dx. 


-/ \yu\^dx<cJ- {\c\'^ + 1) dx + C i J- |VM|"+2d(z: 

dQji(xo) -^Qs^ixo) \"Q 3 r{xo) 

Therefore, the inequality (1!?^ is also satisfied in this case. 

Case 3. Assume Q'^{yo) = 0- 
In this case, inequality (1321) trivially holds. 

In all three cases, the reverse Holder inequality (see Theorem 15.21) shows g G L^ociQ) 
all s G + e) and some e > 0 depending on Rq and n. 



(ii) Higher integrability in the interior. 

This case follows with much less effort and is only sketched here. 

Let xo G H be arbitrary and R > 0 such that Q 2 r{xq) C H. We take a cutoff function 
^ G C^(fl) with 

(a) ^ = 0 in H \ Q 2 R(a;o), (c) C = 1 in Qii(a;o), 

(b) 0<e<linH, (d)|Ve|<|. 

Testing (EH) with ^ = ('^(u — y) and /i = j'Q^jj(xo) yields with the same computation 
as in the case (i): 


f |VMpda;<C f (|c|^ + 1) da; + ^ / |m —/rp dx. 

JQr{xq) JQ2Rixo) ^ JQ2 r{xo) 

The Poincare-Sobolev inequality implies 

n + : 

J- \Vu\^dx<C~[ (|c|'^ + 1) dx + C I -/ iVul'^dxJ 

'^OR(xn) 'J OoTi(xr\^ \^Or>Tf(xn) I 


Applying Theorem 15.21 with g = |Vu|’*+ 2 , q = and / = Cdcj"*^ + l )"+2 finishes the 
proof. ■ 
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6 Existence of weak solutions of ([S^ - logarithmic case 


The challenge here is to establish the integral equation (iii) in Definition l2.3l because the derivative 
of the logarithmic free chemical energy (IA 8 I) becomes singular if one of the c^’s approaches 0. 
We only sketch the proof in this section since all essential ideas can be found in [GarOOl IGarOSb) . 
We use a regularization method suggested in |EL91) and also used in [GarOOl IGarOSb] . 

The energy gradient tensor is assumed to be of the form T = 7 Id with a constant 7 > 0. Define 
a C^(R^) regularization with the regularization parameter (5 > 0 as 

^ 1 
:=0^/(c'=) + -c-Ac, 

fc=l 


with 




X log(x) for d> S, 

X log((5) — I + fy for X < (5. 


Elliott and Luckhaus showed that the regularization is uniformly bounded from below. 

Lemma 6.1 (cf. (EL9l]) There exist constants i5o > 0 and C > 0 such that 

W’='^’^{c)>-C for allcGi:, 5 G {0,So). 


Let qs denote a weak solution in the sense of Definition 12.31 with the free chemical energy = 
p^ch,(5 ^ gy applying Lemma 16.11 and using Gronwall’s inequality in the energy inequality (vi) 
of Definition [231 we can show a-priori estimates analogous as in Section 4 except the a-priori 
estimate of ws ■ 

In the Allen-Gahn case, we have dtcs = —and, consequently, the boundedness of cs in 
and ws £ TT, pointwise lead to boundedness of ws in L^(D;R^). 

In the case of Gahn-Hilliard systems, we can use the following lemma. 

Lemma 6.2 I |Gar00l Lemma 4.3]) There exists a constant C > 0 such that for all 6 £ (0, (5o) 

^ (J VWf’\csit))dx'^ dt<C. 


The proof of this lemma is similar to [GarOOl Lemma 4.3], since all arguments can be adapted 
to our case. Therefore we will omit the proof. 

This lemma and the integral equation 


wsit) dx= [ VWf’\cs{t)) + PIE®i(e(u 5 (t)), csit), zs{t)) dx 
1 dn 


together with the already known boundedness properties shows 















for a constant C > 0. Therefore ws is bounded in T^(0,T; iJ^(O)) by Poincare’s inequality. 
In conclusion, we can extract a subsequence {qs^} such that we have the same convergence 
properties as in Lemma 14.21 As before, we will omit the subscript k. 

Proof of Theorem 12.61 The remaining crucial step is to show that the limit c satisfies Ck > 0 
a.e. on fir for all /c = 1,..., iV and lTc^’^(c 5 ) —>• lT^^^’'°®(c) in as e \ 0. 

To this end, we need an additional boundedness property. 

Lemma 6.3 There exists constants q > 1 and C > 0 such that for all S € (0,(5o) o,i^d all 
k= 

We omit the proof of this lemma, since by utilizing Theorem 15.31 the arguments are analogous to 
[GarOOl Lemma 4.5]. 

Note that 

<5N0 oo otherwise 

holds pointwise a.e. on fir and for all fc = 1,..., iV. Together with Lemma [6.31 we obtain 

> 0 a.e. on TLt 


and 

{(j)^y{cg) —>■ log(c^) + 1 a.e. on 
This and Lemma 16.31 further shows 

^ log(c'') + 1 in L^iflr) 

by Vitali’s convergence theorem. Finally, we can pass to <5 \ 0 in the equation 

f ws - C dxdt = f ^Vcs : VC + PW^’'^(c 5 ) • C + VWy{e{us), cs, zs) ■ C dxdt 
J J six' 

and obtain (iii) from Definition 12.31 

The remaining properties can be easily established as in Section 4. Hence, Theorem l2.6l is proven. 


7 Conclusion 

Materials, which enable the functionality of technical products, change the micro-structure over 
time. Phase separation and coarsening phenomena take place and the complete failure of elec¬ 
tronic devices often results from micro-cracks in solder joints. 

In this work, we have investigated mathematical models describing both phenomena, phase sepa¬ 
ration and damage processes, in a unifying approach. The main aim has been to prove existence 
of weak solutions for elastic Cahn-Hilliard and Allen-Cahn systems coupled with damage phe¬ 
nomena under mild assumptions where the free energy contains 
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• a chemical potential of polynomial or logarithmic type, 

• an inhomogeneous elastic energy, e.g. c, z) = + e)C{c){e — e*(c)) : (e — e*(c)), 

• a quadratic gradient term of the damage variable. 

To this end, several approxmation results have been established as well as different variation¬ 
al techniques, regularization methods and higher integrability results for the strain have been 
applied. 
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